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Abstract 

We report measurements of the cross section and a complete set of polarization 
observables for the Gamow-Teller ^■^C(|?', n)^^N(g.s., 1"*") reaction at a bombarding 
energy of 296 MeV. The data are compared with distorted wave impulse approxi- 
mation calculations employing transition form factors normalized to reproduce the 
observed beta-decay ft value. The cross section is significantly under-predicted by 
the calculations at momentum transfers q > 0.5 fm~^. The discrepancy is partly 
resolved by considering the non-locality of the nuclear mean field. However, the cal- 
culations still under-predict the cross section at large momentum transfers of g ~ 
1.6 fm~^. We also performed calculations employing random phase approximation 
response functions and found that the observed enhancement can be attributed in 
part to pionic correlations in nuclei. 

Key words: Nuclear correlation, Gamow-Teller 1"*" state. Polarization transfer 
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The prediction of pion condensation [1] has prompted extensive experimental 
and theoretical studies of nuclear spin-isospin correlations. Pion condensation 
is expected to occur in cool neutron stars (NS) such as 3C58 [2] because pion 
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condensation can accelerate the cooling of NS [3]. It is believed that pion 
condensation does not occur in normal nuclei; however, precursor phenomena 
may be observed even in normal nuclei if they are in the proximity of the crit- 
ical point of the phase transition. The first proposal for possible evidence of 
a precursor was enhancement of the Ml cross section in proton inelastic scat- 
tering [4,5]. However, this prediction was not supported by measurements of 
^^C(p,p')^^C(l+, T = 1) [6,7]. A possible reason for the absence of the precur- 
sor may be that the Ml cross section involves both pionic (spin-longitudinal) 
and rho-mesonic (spin-transverse) transitions and that the contribution from 
the rho-mesonic transition might mask the pionic effect. 

A further possible source of evidence of a precursor was proposed by Alberico 
et al. [8]. They calculated the pionic and rho-mesonic response functions, Rl 
and Rt, in the quasielastic scattering (QES) region. Their results showed 
significant enhancement in R^/Rt due to nuclear spin-isospin correlations. 
Great effort has been made to extract the spin response functions Rl and 
Rt experimentally in {p,p') scatterings and (p, n) reactions at intermediate 
energies [9]. None of the observed ratios show evidence of the theoretically 
expected enhancement. The fact that the rho-mesonic response Rt is equally 
important in determining the ratio Rl / Rt nieans that the pionic enhancement 
may be masked by the contribution from the rho-mesonic component. Recent 
analysis of QES data shows pionic enhancement in the spin-longitudinal cross 
section that well represents Rl, suggesting that the lack of enhancement in 
Rl/Rt is due to the rho-mesonic component [10]. It should be noted that 
pionic enhancement has been also observed in the pure pionic excitation of 
^^O(p,p')^^O(0-,r = 1) scattering at = 295 MeV [11]. 

Recent progress in the development of high intensity polarized ion sources 
and high efficiency neutron polarimeters has enabled the measurement of a 
complete set of polarization observables for the ^^C(p, n)^^N(g.s., 1+) reaction 
at large momentum transfers covering the critical momentum g ~ 1.7 fm~^ of 
pion condensation. This Gamow-Teller (GT) transition is the isobaric analog 
to the Ml excitation of ^^C(p,p')^^C(l+, T = 1) scattering. In addition, the 
(p, n) reaction is free from isospin mixing effects. Thus, it is very interesting 
to study nuclear correlation effects in this reaction by separating the cross 
section into pionic and rho-mesonic components with polarization observables. 
Furthermore, distorted wave impulse approximation (DWIA) calculations for 
finite nuclei employing nuclear correlations with continuum random phase 
approximation (RPA) are available. 

In this Letter, we present the measurement of the cross section and a com- 
plete set of polarization observables for the ^^C(p, n)^^N(g.s., 1+) reaction. An 
incident beam energy of 296 MeV is used. This is one of the best energies 
to study GT transitions since the spin excitations including GT transitions 
are dominant in the (p, n) reaction near 300 MeV [12]. Furthermore, distor- 
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tion effects become minimum around this incident energy. This allows us to 
extract nuclear structure information reliably by the (p, n) reaction, such as 
the nuclear correlation effects. We compare our results with DWIA calcula- 
tions. Possible evidence of nuclear correlations is observed in the comparison 
between the experimental and theoretical results. We also compare our data 
with DWIA calculations employing RPA response functions including the A 
isobar in order to assess the nuclear correlation effects quantitatively. 

The data were obtained with a neutron time-of-flight (NTOF) system [13] 
with a neutron detector and polarimeter NP0L3 system [14] at the Research 
Center for Nuclear Physics, Osaka University. The NTOF system consists of 
a beam-swinger dipole magnet, a neutron spin-rotation (NSR) magnet, and a 
100-m tunnel. The beam polarization was continuously monitored using two 
p+p scattering polarimeters; its typical magnitude was about 0.70. The beam 
energy was determined to be 296 ± 1 MeV from the kinematic energy shift 
between two peaks from '''Li(p, n)'''Be(g.s. -|- 0.43 MeV) and ^^C(p, n)^^N(g.s.). 
In the beam-swinger system, a beam with a typical current of 500 nA was 
incident on a self-supporting '^'^^C (98.9% ^^C) target with a thickness of 89 
mg/cm^. Neutrons from the target passed through the NSR magnet and were 
measured by the NP0L3 system in the 100-m TOP tunnel with a resolution 
of about 500 keV FWHM. The neutron detection efficiency of NP0L3 was 
determined to be 0.025±0.002 using '^U{p, n)'^Be(g.s. +0.43 MeV) at 0°, whose 
cross section is known for Tp = 80-795 MeV [15]. The neutron polarimetry 
of NP0L3 was calibrated using ^^C{p,ny^N{g.s.) at 0° [14]. The effective 
analyzing power Ay-^s of NP0L3 was determined to be Ay-^s — 0.151±0.007± 
0.004, where the first and second uncertainties are statistical and systematic 
uncertainties, respectively. 

Figure 1 shows the excitation energy spectra of ^^C(]9, n)^^N for momentum 
transfers q = 0.14 fm"^ and 1.7fm-\ The GT 1+ state atE^^O MeV (ground 
state) forms a pronounced peak for q — 0.14fm~^, though it is not fully re- 
solved from the neighboring states for q = 1.7 fm^ . Therefore, we performed 
peak fitting for < 1.5 MeV to extract the yield of the 1+ state. The first 
and second excited states with J'^ = 2+ and 2^ at E^ = 0.96 and 1.19 MeV 
[16] were considered in the peak fitting and were assumed to form a single peak 
because the present energy resolution could not resolve these two peaks. The 
continuum background from wrap-around and ^^C(p,n) events is also consid- 
ered in the peak fitting. The dashed curves in Fig. 1 represent the fits to the 
individual peaks, while the straight dashed line and solid curve represent the 
background and the sum of the peak fitting, respectively. The peak fittings at 
all momentum transfers were satisfactory for extracting the 1"^ yield. 

The differential cross section for ^'^C{p,n)^'^N{g.s.,l^) at Tp = 296 MeV is 
shown in the upper panel of Fig. 2. The data for the analyzing power were also 
measured and the results are displayed in the lower panel. The momentum- 
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transfer dependence was measured in the range q = 0.1-2.2 fm~^, covering the 
maximum at g ^ 1.6 fm^-*^. The data at Tp = 295 McV [17] are also displayed as 
open circles. Both data sets are consistent with each other taking into account 
the statistical and systematic uncertainties. 

We performed DWIA calculations using the computer code CRDW [18]. The op- 
tical model potential (OMP) was deduced from the global OMPs optimized for 
^^C in the proton energy range Tp = 29-1040 McV [19,20]. The free nucleon- 
nuclcon (A^A^) t-matrix parameters were taken from Bugg and Wilkin [21]. The 
single particle wave functions were generated by a Woods-Saxon (WS) poten- 
tial with ro = 1.27 fm-^ oq = 0.67 fm-^ [22], and V^o = 6.5 MeV [23]. The 
depth of the WS potential was adjusted to reproduce the separation energies 
of the 0^3/2 orbits. The transition form factor was normalized to reproduce the 
beta-decay ft value of 13178 s [24], whose corresponding beta-decay strength 
B{GT) was deduced to be B{GT) = 0.873 [25]. The thin solid curves in Fig. 2 
show the DWIA results. The normalization factor for the transition form factor 
is = 0.17. These calculations reproduce the experimental data reasonably 
well for small momentum transfers of g < 0.5 fm~^. However, the angular dis- 
tributions of both the cross section and the analyzing power shift to higher 
momentum transfer. Furthermore, the calculations significantly underestimate 
the cross section at large momentum transfers of g ~ 1.6 fm~^. 

We also investigated the OMP dependence of DWIA calculations by using 
other OMPs [26,27]. The results are shown by the bands in Fig. 2. The OMP 
dependence is clearly seen around the cross section minimum at q ~ 1.4 fm~^. 
However, neither the discrepancy of the angular distributions nor the under- 
estimation of the cross section at large momentum transfer can be resolved. 

We also evaluated the proton-particle and neutron-hole configuration depen- 
dence on the transition form factor. In the above calculations, the Hartrcc- 
Fock (HF) state of ^^C was described as the state fully occupying the Osi/2 
and Op3/2 orbits. Thus, the GT transition is described as a combination of 
a Opi/2 proton-particle and a neutron-hole (0^1/20^3^2)- However, the 

Cohen and Kurath wave functions (CKWFs) [28] show that the 0p3/20p3y^2> 
Op3/20p^/2' ^^^^ 0Pi/20p^/2 configurations also play important roles in the GT 
transition. Therefore we performed DWIA calculations with CKWFs using 
the computer code DWSl [29]. We used the same OMPs and single particle 
wave functions as in the previous calculations. The NN t-matrix parametrized 
by Franey and Love [12] at 325 MeV was used. The dotdash curves in Fig. 2 
denote the DWIA results using the CKWSs with N — 0.89. The cross sec- 
tion is slightly enhanced at all momentum transfers. However, the oscillatory 
structures, especially the momentum transfers where the cross section takes 
a minimum or maximum value, were not improved by using the realistic CK- 
WFs. 
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Therefore, we investigated the non-locahty of the nuclear mean field by using 
a local effective mass approximation [30] in the form of 

m*{r) = mjv - {^^|^| {rriN - m*{0)), (1) 
/ws(Uj 

where rriN is the nucleon mass and fwsir) is a WS radial form. The dashed 
curves in Fig. 2 show the DWIA results with m*(0) = 0.7mAr [31,32] and N 
— 0.17. The angular distributions of both the cross section and the analyzing 
power shift to lower momentum transfer for smaller m*(0), and a standard 
value of m*(0) ~ 0.7mjv improves the agreement with the data. This is because 
the transition form factor effectively moves outwards due to the Perey factor 
effect [33] . However, there is still a discrepancy between the experimental and 
theoretical results around q ~ 1.6 fm~^. 

Next, we consider the nuclear correlation effects. We performed DWIA cal- 
culations with the RPA response functions employing the tt + p + g' model 
interaction Ves- In the spin-isospin channel, Ves can be expressed as [9] 

Kff(g,a;) = Vl(q,a;) + Vr(q,a;) , (2) 

with spin-longitudinal and spin-transverse parts 

VLiq,u;) = Wriq,u;) (n • T2)i<T^ ■ q){cT2 ■ q) 

+ W^\q,u) [{(n ■ T,){a, ■ q){S, ■ q) + {I ^ 2)} + h.c] 
+ Wt\q^ uj) [{(Ti • ■ q)iSl ■ q) + {T, ■ T,){S, ■ q){S, ■ q)} + h.c. 

VT{q,Lj) = W^^{q,u) (ti ■ T2){(Ti x q){(T2 x q) 

+ W^^{q,uj) [{(n ■ T2)(o-i X q){S2 x q) + {1 ^ 2)} + h.c] 

+ W^''{q,Lj) [{(Ti • T^){S, X q){Sl x q) + (Ti • T2)(5i x q){S2 x q)} + h.c. 

(3) 

where Wl and Wt are the spin-longitudinal and spin-transverse strengths, 
respectively, cr and r are the nucleon Pauli spin and isospin matrixes, and S 
and T are the spin and isospin transition operators that excite N to A. The 
strengths Wl and Wt are determined by the pion and rho-meson exchange 
interactions and the Landau-Migdal (LM) interaction Vlm specified by the 
LM parameters g']\[]^, Qna^ 9'aa [9]- For the pion and rho-meson exchange 
interactions, we have used the coupling constants and meson parameters from 
a Bonn potential which treats A explicitly [34]. The LM parameters were 
estimated to be g'j^j^ = 0.65 ± 0.15 and g'j^^ ~ 0.35 ± 0.15 [35] by using the 
peak position of the GT giant resonance and the GT quenching factor sX q — 
[36,37] as well as the isovector spin-longitudinal polarized cross section in 
the QES process at g ~ 1.7 fm^^ [10]. The thick solid curves in Fig. 2 show 
the DWIA result using g'^^^ = 0.65, g'j^^ = 0.35, and m*(0) = 0.7mN with 
N — 0.28. Here, we fixed g'^^ = 0.5 [38] since the g'^^ dependence of the 
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results is very weak. In the continuum RPA, the GT state couples to particle- 
unbound 1+ states, which shifts the response function in the coordinate space 
(r-space) to large r values. Thus the angular distributions further shift to lower 
momentum transfer. The RPA correlation also enhances the cross section at 
large momentum transfers of q ~ 1.6 fm~^. Consequently, the discrepancy 
between the experimental and theoretical results could be resolved in part 
by considering the nuclear correlation effects in the RPA together with the 
non-locality effects of the nuclear mean field. Nuclear correlation effects are 
expected to be more clearly seen in the oscillatory structures of polarized cross 
sections. Therefore, we have separated the cross section into polarized cross 
sections by using the measured polarization observables. 

Figure 3 shows three polarized cross sections, IDq, IDp, and as a func- 
tion of momentum transfer. The spin-scalar polarized cross section IDq is not 
shown because it is very small in the present spin-flip GT case. The error 
bars represent the statistical uncertainties of the data. The systematic un- 
certainties are about 30% for all IDi at g ~ 1.6 fm~^. The solid curves in 
Fig. 3 denote the DWIA calculations with the RPA response function, and 
the bands represent the g'^j^ and g'^^ dependences within g'j^j^ = 0.65=b0.15 
and g'pf^ — 0.35 ± 0.15. The dashed curves are the DWIA calculations with 
the free response function. The momentum transfer dependence of all the 
three spin-dependent IDi shift to lower momentum transfer and are enhanced 
around q = 1.6 fm^^ due to the nuclear correlation effects, which improves the 
agreement with the data. Figure 4 shows the DWIA calculations without the 
normalization. The dashed and dotted curves are the DWIA results with the 
free response function and with the RPA response function employing only 
the LM interaction, respectively. The figure shows quenching of all IDi over 
the whole momentum transfer region due to the repulsive LM interaction, as 
well as forward shifts in the angular distributions due to outward shifts of 
the response functions in r-space by the continuum coupling in the RPA. The 
solid curves are the DWIA results with the RPA response function employing 
the TT + p + g' model interaction Vcs- The pion exchange effects in IDg are 
clearly seen as an enhancement at large momentum transfers. In contrast, the 
rho-meson exchange effects in IDp and ID^ are expected to be small, mainly 
due to the weak momentum transfer dependence of Wt in the present mo- 
mentum transfer region [9]. It should be noted that in the analysis shown in 
Fig. 3, the transition form factors are normalized to reproduce -B(GT), and 
thus the calculated IDi ai'e normalized at small momentum transfers. There- 
fore, the quenching due to the LM interaction is effectively included through 
the normalization factor N at small momentum transfers. However, the mod- 
ification of the angular distributions due to the shape change of the response 
functions in r-space could not be included through this normalization, and 
thus this effect is observed as an enhancement in IDp and IDn in Fig. 3 at 
large momentum transfers oi q ~ 1.5 fm~^. We also note that the momen- 
tum transfer dependence of the NN i-matrix is important for the momentum 
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transfer dependence of IDi. For example, the first minimum of IDq is due to 
the first minimum of the relevant NN t- matrix component [12]. 

The DWIA calculations employing nuclear correlation effects give better de- 
scriptions of all three spin-dependent ID^. However, there are still significant 
discrepancies at large momentum transfers around g ~ 1.6 fm"-*-. The under- 
estimation of the theoretical calculations might be resolved using a realistic HF 
state of ^^C because the DWIA calculation with CKWFs is enhanced around q 
~ 1.6 fm~^ compared with that for the pure Qpi/2Qp^j2 configuration, as shown 
in Fig. 2. Thus, we deduced IDi from the DWIA results with CKWFs; the 
results are shown by the solid curves in Fig. 5. The dashed curves correspond 
to the DWIA results for the pure Qpi/2^Pz}2 configuration. The use of realistic 
CKWFs enhances IDq and IDp at large momentum transfers, and thus the 
discrepancies between the experimental and theoretical results seen in Fig. 3 
would be due to the simplification of the HF state of ^^C in the RPA calcula- 
tions. However, the calculations with CKWFs give smaller IDn values. Thus, 
the reason for the difference in IDn seems to be different to that for IDq and 
IDp. This is possibly due to medium modifications of the effective NN inter- 
action. It should be noted that the enhancement of ID^ at large momentum 
transfers is commonly seen in QES and in stretched state excitations [39] and 
that this enhancement indicates the enhancement of B of the NN scattering 
amplitude in a Kerman-MacManus-Thaler (KMT) representation [40]. Thus, 
the enhancement of B might be responsible for the enhancement of the cor- 
responding rho-mesonic ID^- We note that in both QES and stretched state 
excitations, no modification of the spin-transverse amplitude F has been ob- 
served for /-Dp, whereas a small reduction of the spin- longitudinal amplitude E 
has been indicated for IDq [39]. If we take into account the reduction of -E, the 
enhancement of IDq due to the nuclear correlations should be larger than the 
current prediction using g'^^ = 0.65 and g'^^ = 0.35, which can be achieved 
by using smaller, reasonable LM parameters, g'j^j^ ~ 0.6 and ^r^^ ~ 0.2 [9]. 

In conclusion, our measurement of the cross section and a complete set of 
polarization observables of the ^^C{p,ny^N{g.s.,l~^) reaction enables us to 
study nuclear correlation effects at an intermediate energy of Tp = 296 MeV 
where the theoretical DWIA calculations should be reliable due to the sim- 
ple reaction mechanism. A significant difference in the momentum-transfer 
dependence and an enhancement around q ~ 1.6 fm~^ were observed in the 
cross section compared to the DWIA calculations with Cohen-Kurath wave 
functions. A DWIA calculation employing the RPA response function with 
9'nn — 0-65 and g'j^^ = 0.35 and a local effective mass with m*(0) = O.Tm^v 
reproduces the momentum-transfer dependence and gives an enhancement of 
the cross section at g ^ 1.6 fm~^. All three spin-dependent polarized cross 
sections also support the existence of nuclear correlation effects. In theoreti- 
cal calculations, the use of a realistic shell-model wave function seems to be 
important to reproduce the experimental data at large momentum transfers. 
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However, this effect could not be included in the present RPA calculations. 
Therefore, a more comprehensive and detailed theoretical analysis is needed 
to explain quantitatively the nuclear correlations inside nuclei. 
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Fig. 1. Excitation energy spectra for ^^Cip,n)^^N at Tp = 296 MeV and q = 
0.14fm~^ (upper panel) and q = 1.7fm~^ (lower panel). The dashed curves and 
straight dashed line represent fits to the individual peaks and background, respec- 
tively. The solid curve shows the sum of the peak fitting. 
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Fig. 2. Comparison between experimental and theoretical cross sections (upper 
panel) and analyzing powers (lower panel) for ^^C(p, n)^^N(g.s., 1+) at Tp = 296 
MeV. The thin solid curves are the DWIA calculations for a 0pi/20p^/^2 configuration 
and the bands represent the OMP dependence. The dotdash curves denote the 
DWIA results with CKWFs. The dashed and thick solid curves are, respectively, 
the DWIA results with the free response function employing m*(0) = 0.7m at and 
with the RPA response function employing g'^j^ = 0.65, g'j^^ = 0.35, and m*(0) = 
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Fig. 3. Comparison between experimental and theoretical IDq (left panel), I Dp 
(middle panel), and ID^ (right panel) for ^'^C{p,nY'^^{g.s., 1+) at Tp = 296 MeV. 
The solid curves are the DWIA calculations with the RPA response functions using 
(Snn^ 9na^ 9'aa) ~ (0.65,0.35,0.50) and m*(0) = 0.7mN- The bands represent the 
g'ji^jif and g'ji^/^ dependences of the DWIA results within 0.50 < (7jviV — 0-80 and 
0.20 < 

9'na — 0.50. The dashed curves show the DWIA results with the free response 
functions using m*(0) = 0.7mN. 
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Fig. 4. Theoretical calculations without normalization for IDq (left panel), I Dp 
(middle panel), and (right panel) for '^'^C{p,nY'^^{g.s.,l+) at Tp = 296 MeV. 
The dashed curves show the DWIA results with the free response function. The 
dotted and solid curves are the DWIA calculations with the RPA response func- 
tions employing the Landau-Migdal interaction [(f) only and the it + p + c/ model 
interaction, respectively. 
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Fig. 5. DWIA predictions for '^'^C{p,nY'^^{g.s.,l+) at Tp = 296 MeV. The solid 
curves are the DWIA results with CKWFs. The dashed curves represent the DWIA 
results for a ^Pi/2^P^j2 configuration. The experimental data are same as those in 
Fig. 3. 
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